SUBMITTED TO ApJ 

Preprint typeset using LAT^X style emulateapj v. 14/09/00 



THE INTERNAL STRUCTURE AND FORMATION OF EARLY-TYPE GALAXIES: 
GRAVITATIONAL-LENS SYSTEM MG2016+112 AT Z= 1.004 



THE 



(N 
O 
O 
(N 

u 

< 
O 



(N 
> 
(N 

co 

(N 
O 
(N 
O 

Or 

6 



TOMMASO TREU 

California Institute of Technology, Astronomy, mailcode 105-24, Pasadena, CA 91125 

Leon V.E. Koopmans 

California Institute of Technology, Theoretical Astrophysics Including Relativity (TAPIR), mailcode 130-33, 

Pasadena, CA 91125 
submitted to ApJ 

ABSTRACT 

We combine our recent measurements of the velocity dispersion and the surface brightness profile of 
the lens galaxy D in the system MG2016+112 (z = 1.004) with constraints from gravitational lensing to 
study its internal mass distribution. We find that: (i) dark matter accounts for >50% of the total mass 
within the Einstein radius (99% CL), whereas ^75% is the more likely contribution. In particular, we 
can exclude at the 8-cr level that mass follows light inside the Einstein radius with a constant mass-to- 
light ratio (M/L). (ii) the total mass distribution inside the Einstein radius is well-described by a density 
profile oc r -7 with an effective slope 7' = 2.0 ±0.1 ±0.1, including random and systematic uncertainties, 
(iii) The offset of galaxy D from the local Fundamental Plane independently constrains the stellar M/L, 
and matches the range derived from our models, leading to a more stringent lower limit of >60% on the 
fraction of dark matter within the Einstein radius (99% CL). 

Under the assumption of adiabatic contraction, we show that the inner slope of the dark matter halo 
before the baryons collapsed to form the lens galaxy is % < 1.4 (68% CL), only marginally consistent with 
the highest-resolution cold dark matter simulations that indicate 7i~1.5. This might indicate that either 
adiabatic contraction is a poor description of early-type galaxy formation or that additional processes 
play a role as well. Indeed, the apparently isothermal density distribution inside the Einstein radius, is 
not a natural outcome of adiabatic contraction models, where it appears to be a mere coincidence. By 
contrast, we argue that isothermality might be the result of a stronger coupling between luminous and 
dark-matter, possibly the result of (incomplete) violent relaxation processes during the formation of the 
innermost regions of the galaxy. Hence, we conclude that galaxy D appears already relaxed ~8 Gyr ago. 
We briefly discuss the importance of our results for lens statistics and the determination of the Hubble 
Constant from gravitational-lens time delays. 
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1. INTRODUCTION 

Understanding the internal structure of galaxies, is cru- 
cial to explain the variety of their observed morphologies 
and the physical processes that created them. It also pro- 
vides an opportunity to test physics on galactic scales. 

Because early-type (E/S0) galaxies are relatively simple 
systems, dynamical modeling of their internal structure 
can be attempted with hope for success (e.g. de Zeeuw & 
Franx 1991; Bertin & Stiavelli 1993; Merritt 1999). Knowl- 
edge of this internal structure provides a vital discriminant 
to the processes involved in their formation (e.g. van Al- 
bada 1982), which is currently a highly controversial issue 
(e.g. Schade et al. 1999). 

In particular, in the cold dark matter (CDM) cosmolog- 
ical model, galaxies are embedded in dark matter halos, 
which dominate the total mass of the system and have a 
characteristic mass distribution (Navarro, Frenk & White 
1997, hereafter NFW; Moore et al. 1998; Bullock et al. 
2001; Ghigna et al. 2000). Measuring the distribution of 
dark matter in E/S0 galaxies would therefore not only shed 
light on their formation history, but also provide a pow- 



erful test of the CDM cosmological model. This test is 
similar, but complementary in the scales and physical con- 
ditions probed, to the still controversial tests performed 
using the kinematics of low surface-brightness and dwarf 
galaxies (e.g. McGaugh & de Blok 1998; van den Bosch 
et al. 2000; Swaters, Madore & Trewhella 2000; Salucci & 
Burkert 2000; Borriello & Salucci 2001; van den Bosch & 
Swaters 2001; de Blok et al. 2001; de Blok & Bosma 2002; 
Jimenez, Verde & Peng 2002). 

Unfortunately, E/S0 galaxies generally lack simple dy- 
namical tracers at large radii, such as HI in spiral galaxies, 
and little is known about their kinematics beyond the re- 
gion dominated by stellar mass (1-2 effective radii, R e ). 
Hence, the amount and distribution of dark matter within 
E/S0 galaxies remain to date poorly constrained even in 
the local Universe (e. g. Bertin et al. 1994; Franx, van 
Gorkom & de Zeeuw 1994; Carollo et al. 1995; Rix et al. 
1997; Gerhard et al. 2001). In a few exceptional cases 
kinematic tracers at large radii have been found and used 
to investigate the presence and distribution of dark matter 
halos (e. s. Mould et al. 1990; Franx et al. 1994; Hui et 



Based on observations collected at W. M. Keck Observatory, which is operated jointly by the California Institute of Technology and the 
University of California, and with the NASA/ESA Hubble Space Telescope, obtained at STScI, which is operated by AURA, under NASA 
contract NAS5-26555. 
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In a previous paper (Koopmans & Treu 2002; hereafter 
KT02) we described the Keck and HST observations of 
the E/SO lens galaxy D in MG2016+112, the most dis- 
tant spectroscopically confirmed lens galaxy known (z = 
1.004). A review of the current observations of this system 
can be found in Koopmans ct al. (2002; hereafter K02), 
who also propose a detailed lens model, where the im- 
age separation is caused predominantly by galaxy D, and 
find a central velocity dispersion <r=320-340 kins" 1 for 
an isothermal total mass distribution. This value should 
correspond closely to the central stellar velocity dispersion 
(e. g. Kochanek 1994; Kochanek et al. 2000). Our direct 
measurement (a = 328 ± 32 kms -1 ; KT02) shows that D 
is indeed a massive old E/S0 galaxy, possibly embedded 
in a proto-cluster, with a central stellar velocity dispersion 
in agreement with that inferred from the lens models. 

In this paper we present a detailed analysis of the inter- 
nal structure of galaxy D, based on the observed velocity 
dispersion and the exquisite HST images available from 
the CASTLES database. In particular, we focus on the 
density profiles of (i) the dark matter halo and (ii) the 
total luminous plus dark mass distribution. Apart for its 
high rcdshift (corresponding to ^ 8 Gyr look-back time 3 ), 
this system is particularly interesting because its Einstein 
radius (R,Einst=13.7±0.1 kpc) is much larger than the ef- 
fective radius of the lens galaxy (R c = 2.7 ± 0.5 kpc). For 
this reason, we can use gravitational lensing to probe the 
mass distribution to much larger radii than is typically 
possible, even in local E/S0 galaxies. 

The paper is organized as follows. First, the lumi- 
nous and dark-matter models are described in Sect. 2. In 
Sect. 3, constraints are placed on the dark and luminous 
matter distributions inside the Einstein radius, and com- 
pared to CDM predictions in Sect. 4. In Sect. 5, we use the 
results found for MG2016+112 to predict observables for 
other similar systems. Conclusions and a discussion are 
given in Sect. 6. Throughout this paper, r is the radial 
coordinate in 3-D space, while R is the radial coordinate 
in 2-D projected space. 

2. MODELS 

Based on the lens models in K02, one finds a total mass 
enclosed within the Einstein radius, M(< i?Einst) = Me = 
1.1 x 10 12 M Q . The enclosed mass has been corrected for 
the slight cllipticity of the lens potential and, because the 
source lies almost on the observer-lens axis, is nearly in- 
dependent of the mass distribution assumed for the lens 
model (e.g. Schneider, Ehlers & Falco 1992; Chapter 8). 
An error of ~ 10% is assumed for the enclosed mass (K02). 
This error is larger than what one might expect from a 
range of lens models (e.g. see Kochanek 1991), but is cho- 
sen to include any remaining uncertainties related to pos- 
sibly undetected mass contributions from nearby galaxies. 
Note, however, that to first order, the contribution from 
the most massive nearby structures has been corrected for 
in the lens models (K02). The error estimate on the en- 
closed mass can thus be regarded as conservative. Later, 
we will estimate the effect of this potential systematic error 
on our conclusions and show it to be negligible. 

2 see 'www.astro.caltech.edu/~tt/LSD' or 'www.its.caltech.edu/~koopmans/LSD' 

3 throughout the paper, we assume the Hubble Constant, the matter density, and the cosmological constant to be Ho=65 kms -1 Mpc -1 , 
f! m = 0.3 and Q\ = 0.7, respectively. 



al. 1995; Arnaboldi et al. 1996). However, these methods 
do not seem likely to be applicable at significant look-back 
times and thus to the investigation of the evolution of the 
internal structure of E/SO galaxies. 

An independent measure of the mass distribution of 
E/SO galaxies is provided by gravitational lensing, which 
is not sensitive to its dynamical state. For example, gravi- 
tational lensing alone has been used to show the existence 
of dark matter in individual E/SO galaxies (e. g. Kochanek 
1995) and recently to constrain the dark matter halo pro- 
file from gravitational-lens statistics (e.g. Keeton 2001; 
Kochanek & White 2001), and galaxy-galaxy lensing sig- 
nal (e.g. Seljak 2002). The combination of dynamical and 
lensing analyses can significantly reduce the degeneracies 
inherent to each method and result in tighter constraints 
on the dark-matter distribution, as well as on the luminous 
and total mass distribution. In addition, most lenses are 
found at intermediate redshift (0.1 < z < 1) and there- 
fore offer the possibility to probe the mass distribution at 
significant look-back times. 

1.1. The Lenses Structure and Dynamics (LSD) Survey 

To investigate the internal structure of E/SO galaxies 
and test CDM predictions, we have started a program to 
measure the kinematic profiles of a sample of E/SO galaxies 
that are also gravitational lenses; the Lenses Structure and 
Dynamics Survey (hereafter the LSD Survey 2 ). By means 
of deep spectroscopy, using the Echelle Spectrograph and 
Imager (ESI) at the W. M. Keck-II Telescope, we are de- 
riving extended kinematic profiles (beyond the effective 
radius) for the nearest and brightest lenses, along major 
and sometimes minor axes, while for the most distant and 
faint ones we are measuring a luminosity-weighted veloc- 
ity dispersion. The ultimate aim of the LSD Survey is to 
produce a database of measurements of internal kinemat- 
ics for the largest possible number of E/SO lens galaxies, 
covering a range of redshifts and masses. 

The target lenses have to: (i) have E/SO morphology, as 
determined from Hubble Space Telescope (HST) images; 
(ii) be relatively undisturbed by the light of the lensed ob- 
ject at ground based resolution; (iii) be relatively isolated, 
since a nearby cluster, large group, or massive companion, 
would introduce additional uncertainties in the measure- 
ment of the mass of the lens; (iv) be bright enough that 
accurate internal kinematics can be determined in a few 
hours of integration. A preference is given to objects where 
the redshift of the background source is known. If the red- 
shift of the source is not known and the source is bright 
enough that the redshift can be determined in a few hours 
of integration time, the redshift is also determined. 

In addition to help clarifying the origin of E/SO galax- 
ies, other important results will be obtained if a better 
understanding of the internal structure of E/SO and its 
evolution with rcdshift is achieved. For example we can 
mention the determination of the Hubble constant from 
time-delay, and the calculations of lens cross-sections. 

1.2. A case study: MG2016+112 



The internal structure of MG2016+112 
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The total luminosity of the system is (1.6±0.2)10 11 L Qi s 
(KT02), implying an average mass-to- light ratio M/Ls = 
8 ± 1 Mq/Lq.b inside the Einstein radius, assuming an 
i? 1 / 4 surface-brightness profile. This value is similar to 
the values found for the stellar populations of local E/SO 
galaxies (e.g. Gerhard et al. 2001). However, considering 
the much smaller age of the stellar populations of galaxy 
D - estimated to be 3.0 ± 0.8 Gyr and not larger than the 
age of the Universe at z = 1.004 (~ 6 Gyr; KT02) - the 
stellar M/L ratio of galaxy D must be smaller than in the 
local Universe, due to luminosity evolution of the stellar 
population. Given the large value of M/L, we conclude 
that a significant amount of dark matter is likely to be 
present inside i?Einst- 

We choose to describe the mass distribution of 
MG2016+112 with two spherical components, one for the 
luminous matter and one for the dark matter halo. In par- 
ticular, we model the luminous mass distribution with a 
Hernquist (1990) model 



M*r« 



2irr(r + r») 3 ' 



(1) 



that well reproduces the i? 1 / 4 surface brightness profile for 
r* = i? e /1.8153 (M* is the total stellar mass) and is sim- 
ple to treat analytically (e.g. Ciotti, Lanzoni & Renzini 
1996). The dark-matter distribution is modeled as 



Pd(r) = 



Pd,0 



(r/r b )t{l + {r/r b yy 



(2) 



which closely describes a NFW profile for 7 = 1, and has 
the typical asymptotic behavior at large radii found from 
numerical simulations of dark matter halos oc r~ 3 (e. g. 
Ghigna et al. 2000). Note that, for r b significantly exceed- 
ing i?Einst, the particular choice of the outer slope has only 
a negligible effect on the kinematics of the lens galaxy. In 
particular, for r& — > 00 and 7 = 2, Eq. 2 reduces to a sin- 
gular isothermal sphere, a functional form often used to 
describe dark matter halos. In addition, Eq. 2 provides 
a simple parametrization of the controversial inner slope 
for r < Tb, the region better constrained by our observa- 
tions. The length scale (rb) and the density scale (pd,o) 
determine the virial mass of the dark matter halo (e. g. 
Bullock et al. 2001). 

The velocity dispersion is computed solving the spheri- 
cal Jeans equation (e.g. Binney & Tremaine 1987), assum- 
ing an Osipkov-Merritt (Osipkov 1979; Merritt 1985a,b) 
parametrization of the anisotropy (3: 



al r 2 



(3) 



3. LUMINOUS AND DARK MATTER IN MG2016 + H2 

3.1. Constraints from kinematics, tensing geometry, and 
surface brightness distribution 

The effective radius and the total mass (Me) within 
the Einstein radius are fixed by the observations (KT02). 
This leaves four free parameters in the model: the in- 
ner slope (7) of the dark matter halo, the length scale of 
the dark matter component (rb), the mass-to- light ratio 
of the luminous component (WL^/Lb) and the anisotropy 
radius (ri). In Fig. 1 we show the likelihood contours 4 
of 7 and M*/Lb, given the observed velocity dispersion 
{o~ap = 304 ± 27 krns" 1 ) for a range of values of rb- Since 
changing r, has very little effect on the likelihood contours, 
only the ones obtained for = R e are shown (for more dis- 
cussions see also also the rest of this section and in Sect. 5). 



r t /R EI „, t =l 



i-b/R E i„„=2 




where o~e and ay are the tangential and radial component 
of the velocity dispersion. Note that (3 > by definition 
and therefore tangentially anisotropic models are not con- 
sidered in this analysis. Finally, the velocity dispersion 
as observed through our spectroscopic aperture (= a ap ; 
effectively corresponding to a circular aperture of radius 
0'.'65 ~ 2R e ) is computed as the average of the projected 
velocity dispersion weighted by the surface brightness. 

4 Gaussian error distributions are assumed. 

5 We use the term effective slope to emphasize that this is the slope that we get with a single power law. The real mass distribution could 
change slope between the constraints given within the spectroscopic aperture radius and at REinst, as l° n g as the effective trend is preserved. 

6 defined as the luminosity- weighted velocity dispersion within a circular aperture of radius R c /8 



Fig. 1. — Constraints on the inner slope of the dark matter halo 
(7) and the stellar mass-to-light ratio (M, / Lb , in solar units) . Thin 
contours represent the 68%, 95%, 99% confidence limits given Mb, 
R e and cr ap . The thick contours represent the same limits, if the 
independent measurement of M t /Lb from the evolution of the FP 
is included (Eq. 6). The subpanels present the results for different 
values of the dark-matter length scale (r;,). The contours change 
negligibly for different values of the anisotropy radius anisotropy 
radius (ri/R e = 0.5,1,2), thus only the case ri/R e = 1 is shown. 
See Sect. 3 for a detailed discussion. 

Two main conclusions can be drawn from Fig. 1. First, 

since M*/L B < 4M /L aB (99% CL), i.e. smaller than 
the total M/Lb~S ± IMq/Lq^, dark matter is present 
within i?Einst and contributes > 50% (99% CL) of the 
mass. As an independent check, if we assume that Me 
arises from the luminous component only with a constant 
M^/Lb = 8M0/L0 ; s, we find a ap — 530 kms -1 (virtually 
independent of the isotropy radius), inconsistent with the 
observed value at better than 8a. Second, we find that 
for large values of r & /i? E inst (e. g. r b /R Einst = 2 shown in 
Fig. 1) and small values of M^/Lb, 7 approaches the value 
2. Indeed, results from calculations with M*/Lb — > and 
r b — ► 00, show that 7 = 2.0 ± 0.1 (1-a error). In other 
words, the total mass can be very well described by a 
single power-law density distribution p t oc r -7 with ef- 
fective 5 slope 7' = 2.0 ± 0.1. A systematic error of 10% 
on Me changes the inferred slope by only 0.05, whereas 
a 20% error (1-ct) on R c affects the result by only 0.02. 
Similarly, a 20% error (1-ct) on the assumed aperture ra- 
dius of 0'.'65 changes the result by 0.03. Even a factor 10 
deviation from j-j = R c results in < 0.05 error, although 
the effect on the central velocity dispersion 6 a in that case 
is much more profound, with significantly larger values of 
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a for smaller rj (i.e. more radial orbits, see Fig. 4). Hence, 
7' is robust against reasonable changes in the parameters 
rj, i? c , Me, the aperture radius, and r& > i?Einst- 

We can make a very simple argument to understand why 
the effective slope is so tightly constrained by our obser- 
vations. For a total mass profile py oc r~ 7 the enclosed 

mass is M(< r) oc r 3 1 . Therefore, the mass within the 
effective radius (M e ) depends on the mass within the Ein- 
stein radius (Me) as 



M e = M E 



Re 



R 



Litis!. 



3-7' 



= M EX y\ 



(4) 



where xe = REinst/Re has been introduced for simplicity 
of notation. Hence, assuming M e oc a 2 for a given R e , 



5i 



25a 



(5) 



a\\og{x E )Y 

where 5j' and 5a are the errors on the effective slope 
and luminosity-weighted velocity dispersion, respectively. 
Using the values for MG2016+112 in Eq. 5 we recover 
5-y' « 0.1 in good agreement with the error inferred from 
our model. 

3.2. Further constraints: the evolution of the 
Fundamental Plane 

Recent studies have shown that E/SO galaxies both in 
clusters and in the field define a tight Fundamental Plane 
(FP; Dressier et al. 1987; Djorgovski & Davis 1987) out to 
z ~ 0.7 - 0.8 (van Dokkum et al. 1998; Treu et al. 1999, 
2002; see also Kochanek et al. 2000), with slopes very sim- 
ilar to the ones observed in the local Universe. Assuming 
that galaxy D in MG2016+112 lies on a FP with slopes 
as in the local Universe, we can obtain an additional con- 
straint on its internal structure. In fact, the evolution of 
the intercept of the FP with rcdshift can be related to 
the evolution of the average effective mass to light ratio 
(A log M/Lb] see e. g. Treu et al. 2001 for discussion). 
For D we obtained A log M/Lb = -0.62 ±0.08, i.e. inter- 
mediate between the cluster and field value (KT02). We 
can use this measurement to infer M*/Lb of D assuming 
that 

log(M,/L B ) z - Iog(M,/L B )o + A log(M/L s ), (6) 

where the second term on the right hand side of the 
equation is measured from the evolution of the FP. The 
first term on the right hand side of the equation can be 
measured for local E/S0 galaxies: for example using the 
data from Gerhard et al. (2001; and references therein), 
we find a good correlation between a 2 R e and M^/Lb- 
Using the values of MG2016+112, this correlation yields 
M*/L B {z = 0) = 7.3 ± 2.1M /L 0) b, where the error in- 
cludes the contribution from the scatter between dynam- 
ical and stellar population measurements 7 . Hence, from 
Eq. 6 we find M*/L B {z = 1.004) = 1.8 ± O.7M /L 0jS . 
We note two things: (i) the value M*/Lb = 8M /L 0i b 
required to explain the enclosed mass solely by luminous 
matter is inconsistent with the above, independently- 
derived, value at the 9-a level, and (ii) the allowed 
range of M*/Lb found from our dynamical modeling 
(Fig. 1) overlaps very nicely with the value derived from 
the evolution of the FP. Using this additional constraint 



we obtain the probability contours plotted as thick solid 
lines in Fig. 1. Again, the results depend very mod- 
estly on the anisotropy radius, while the slope 7 be- 
comes steeper as r& increases. In general, 7 = 2 can 
be ruled out at better than 2-4 a (depending on r^), 
and therefore the dark matter profile is not isothermal. 
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Fig. 2. — Dark and luminous matter best-fit density (in units of 
10 10 Mq kpc -3 ) and velocity dispersion profiles (in kms -1 ) as func- 
tion of radius (in kpc). Lower panel: Density profile of the luminous 
component (dot-dashed), of the dark matter halo (dotted) and the 
sum of the two (solid) for several characteristic best-fitting models. 
An r~ 2 density profile (dashed) reproduces the total density profile 
remarkably accurately. Middle panel: Relative residuals of the to- 
tal mass distribution (ptot) with respect to power laws py = r — 7 , 
with 7' = 2.2, 2, 1.8 from top to bottom. Upper panel: Luminosity- 
weighted velocity dispersion profile as a function of circular aper- 
ture. The solid lines are obtained for the best-fitting two-component 
model. The thin dot-dashed lines correspond to power-law mass 
distributions with 7' =2.2,2,1.8 from top to bottom. The error bar 
represents our measurement a ap = 304 ± 27km s _1 . We assume 
M*/L B = 1.8M /L 0iS . 



Let us now examine the total mass distribution of the 
two-component model in more detail with the help of the 
quantities plotted in Fig. 2. First (lower panel), we show 
the density profiles of a characteristic best-fitting model 
{r b = 2R e , n = R e ; M*/L = 1.8 M /L 0)S ). Even though 



luminous and dark matter themselves are not isothcrmally 

7 A very similar result is obtained by taking simply the average stellar M/L of the sample 7.8 ± 2.7M0/Lq i b. 
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distributed, the total mass distribution is very close to 
isothermal, confirming what we found in Section 3.1. Sec- 
ond (middle panel), we show the relative difference be- 
tween a power law mass profile py oc r 1 and the total 
mass profile obtained with our two-component model. A 
profile with 7' = 2 describes the total mass profile very 
well. By contrast, 7' = 1.8 or 7' = 2.2 result in sig- 
nificantly worst representations of the total mass profile. 
Third (upper panel), we compare the luminosity- weighted 
velocity dispersion profile of our two-component model 
with those obtained for 7' =1.8, 2.0 and 2.2. The 7' = 2 
power law mass profile and the two-component model pro- 
duce almost identical velocity dispersion profiles, while 
7' = 1.8, 2.2 are cleary inconsistent with the observed a ap 
(point with error bar in the plot). For completeness, we 
note that we find cr/a ap = 1.06, consistent with our as- 
sumption of 1.08 ± 0.05 in KT02. 

Finally, we check that our best-fitting models corre- 
spond indeed to physical solutions in the sense that they 
can be generated by a positive distribution function (i.e. 
they are consistent models, see e.g. Ciotti 1999 and refer- 
ences therein for discussion). For our choice of values of 
Ti/R ei a luminous Hernquist model embedded in a total 
isothermal mass distribution always satisfies the sufficient 
condition given in Eq. 6 of Ciotti (1999), and therefore our 
best-fitting models are consistent. 

4. COMPARISON WITH CDM MODELS 

In order to make a proper comparison with CDM pre- 
dictions, we have to consider that the universal profiles 
of dark matter halos can be significantly modified by the 
formation of a galaxy within the halo. A simple treat- 
ment of this phenomenon can be given by assuming the so- 
called "adiabatic approximation" (Blumcnthal et al. 1986; 
see also Mo, Mao & White 1998). Under this approxi- 
mation, for any given initial mass distribution (baryons 
are assumed to follow dark matter at the beginning) we 
can calculate the final dark matter distribution after the 
galaxy has collapsed, for example to an Hernquist mass 
distribution (e. g. Keeton 2001; see also Rix et al. 1997). 
Here, we use the approximation backwards to recover the 
initial dark matter density profile. In order to facilitate 
comparison to CDM predictions, we fit the same functional 
form as in Eq. 2 to the initial density profile by minimizing 
the fractional difference in density integrated over 0-100 
kpc, and derive the initial length scale Tb.i and an initial 
slope 7, . Although it was not guaranteed a priori, we find 
that this form gives a reasonable description of the initial 
density profile as well. 

Before proceeding to the calculation, we reduce the de- 
grees of freedom of the problem by noticing that CDM 
models indicate that the length scale of systems as massive 
as galaxy D should be significantly larger than ~ 15 kpc, 
even at z ~ 1 (Bullock et al. 2001). Therefore, in the spirit 
of the comparison with CDM predictions, we can limit our 
parameter space to r& > i?Einst> where the results depend 
very little on the exact value of rf,. The length scale does 
not change significantly (< 30% depending on the slope 
for ru < 40 kpc) during collapse, therefore our assump- 
tion on r;, being large because r^i is large is justified in 
this context. For definiteness, we pick r\, = 2i?Einst so 
that the observed profiles are shown in Fig. 1 and 2. With 



this assumption, the initial slope 7$ depends only on the 
observed slope 7 and on M*/Lb- 

In Fig 3 we show the contour levels of 7, in the 
7 - M*/L B plane. As M*/L B —> , the effect of 
baryons becomes negligible and 7 — ► 7$. On the other 
hand, as M^/Lb increases, the collapse becomes in- 
creasingly important and flatter initial dark matter ha- 
los correspond to increasingly steeper final halos. To 
constrain 73 with the observations, we overlay the likeli- 
hood contours of 7 — M^/Lg derived from our dynamical 
model (Fig. 1), and consider only the portion of the 7, 
contours enclosed by the likelihood contours. Consid- 
ering only monotonically decreasing profiles, the range 
7i = 0.0 — 1.4 maps onto a very tight range of final slopes 
7 « 1.30 — 1.75. In conclusion - within the context 
of this simple adiabatic contraction model -74 = 1.5 
is only marginally consistent with our observations, 
whereas ji = 1 fits within the 68% likelihood contour. 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 
M./L B 



Fig. 3. — Relation between the final (7) and the initial (7,) inner 
slope of the CDM dark-matter halo of galaxy D in MG2016+112, 
and the stellar mass-to-light ratio (M„/Lb), based on a simple adi- 
abatic contraction model (see text). The loci corresponding to dif- 
ferent initial slopes are shown, labeled by 7$. In particular, the 
two thick contours labeled 1.0 and 1.5 (=7i) indicate the NFW 
and Moore et al. (1998) profiles, respectively. The probability 
contours given the observations of MG2016+112 for a model with 
= 2 i?Einst and ri = R c (sec Fig.l) arc ovcrplotted and labeled by 
the probability level. As in Fig. 1 solid lines include the constraint 
from the FP evolution, dashed lines do not. 

5. KINEMATICS OF HERNQUIST LUMINOSITY DENSITY 
PROFILES IN A LOGARITHMIC GRAVITATIONAL 
POTENTIAL 

In Sect. 3 we have shown that the total mass distribu- 
tion of galaxy D is well described by an isothermal sphere 
in the range ~ 1-15 kpc. Here, we generalize the results 
found for galaxy D to compute the velocity dispersion pro- 
files that we expect to find for other lens E/S0 galaxies. 

Let us assume that E/S0 galaxies have total mass pro- 
files that follow ptot °c r~ 2 between a fraction of an 
effective radius and a few effective radii. Solving the 
Jeans Equation for a Hernquist luminosity-density pro- 
file embedded in a logarithmic potential, we can com- 
pute the velocity dispersion profile, depending only on the 
effective radius (R e ), the "Einstein velocity dispersion" 
(oEinst = \/G M E /R-Einst), and the anisotropy radius (r») 
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In Fig. 4, wc show the luminosity-weighted velocity 
dispersion profile (in units of R e and OEinst), f° r sev- 
eral values of r^. Note that the profiles approach the 
same curve for large R/R e independently of j-j. Hence, 
the effect of changing on the luminosity-weighted ve- 
locity dispersion beyond the effective radius, is not as 
dramatic as in the inner portions of E/SO galaxies. For 
this very reason, the luminosity-weighted velocity disper- 
sion does not firmly constrain r,. On the other hand, 
the velocity dispersion profile itself is much more sensi- 
tive to Ti. Therefore, a spatially resolved measurement of 
er(r), extending from a fraction of to a few R e , will allow 
us to constrain both n and the total mass distribution. 
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Fig. 4. — Luminosity-weighted velocity dispersion profile for a 
Hcrnquist luminosity-density profile in logarithmic potential. The 
axes are scaled in units of R e and CTEinst (see text). From top to 
bottom curves are shown for ri/R e =0.2,0.5,l,2,oo (where oo cor- 
responds to isotropic models). Note the smaller effect of changing 
anisotropy radius r;, as compared to changing the effective slope 7' 
of the total mass distribution (upper panels in Fig. 2; to facilitate 
comparison the two figures show the same range of velocity disper- 
sion and radius). Our measurement of MG2016+112 is shown as a 
filled circle. The value for 0047-0281 using data from Warren et al. 
(1998) and Kochanek et al. (2000) is shown as an open triangle. 

To our knowledge, the only system (besides 
MG2016+112) satisfying the LSD requirements 8 (Sec 1.1), 
for which velocity dispersion of the lens galaxy has been 
measured is 0047-281 (Warren et al. 1996, 1998; the 
effective radius and image separation are measured by 
Kochanek et al. 2000). As shown in Figure 4, the ob- 
served velocity dispersion of the lens galaxy in 0047-281 
agrees very well with the value computed with our model. 
In future papers we will present velocity dispersion (pro- 
files) for other systems, and further test this model. 

6. SUMMARY & DISCUSSION 

6.1. Summary 

We have shown that the combination of constraints from 
stellar kinematics and gravitational lensing is extremely 



powerful in determining the mass profiles of luminous and 
dark matter in E/S0 galaxies even at large cosmological 
distances. In particular, for galaxy D in MG2016+112 
(z = 1.004) we have demonstrated that: 

(i) Dark matter is unambiguously present and accounts for 
>50% (99% CL) of the mass within the Einstein radius. A 
mass distribution following light with constant M/L can 
be ruled out at more than 8-cr level. Considering an addi- 
tional constraint on M*/Lb from the evolution of the FP, 
the smallest allowable fraction of dark matter increases to 
60% (99% CL). 

(ii) The total mass distribution (dark+luminous) inside 
the Einstein radius is well described by an isothermal 
sphere, i.e. p tot oc r~ 7 with 7' = 2.0 ± 0.1. This result 
is robust and does not depend significantly on the dynam- 
ical state of the system or on the modeling of the mass 
distribution (one or two component models). The result 
is only weakly dependent on errors on Me, ri, R e and the 
aperture radius, which in quadrature add at most another 
0.1 to the error. 

(iii) Using two-component (dark and luminous) spherical 
models, we find that the inner slope of the dark matter 
halo 7 < 2.0 at the 95% CL, i. e. the halo is not isother- 
mal inside the Einstein radius. 

(iv) Using a simple adiabatic contraction model, we re- 
late the observed inner slope 7 to the inner slope ji of 
the dark matter halo before collapse. In particular we find 
the upper limit 7; < 1.4(2.0) at 68(99)% CL, marginally 
consistent with the high resolution CDM simulations that 
indicate 7i~1.5 (Ghigna et al. 2000; but see Power et al. 
2002 for a discussion of the related uncertainties). 

These results are robust because the stellar velocity dis- 
persion and gravitational lensing provide two independent 
mass measurements at two well-separated radii, the effec- 
tive radius (2.7 kpc) and the Einstein radius (13.7 kpc 
~ 5R C ). This leads to a very well-defined slope inside the 
Einstein radius, and to stringent limits on the fraction of 
dark matter. 

Higher signal to noise - and possibly spatially resolved 
observations - are needed to improve these constraints. 
However, this is unlikely to be feasible with the current 
generation of space or ground-based telescopes and will 
probably need the Next Generation Space Telescope or 
large ground-based telescopes with adaptive optics. Nev- 
ertheless, when the measurement and analysis of the dozen 
systems targeted by the LSD Survey will be completed, wc 
expect to be able to obtain more stringent constraints by 
combining the results from the individual systems. 

Finally, in Sect. 5 we have provided a simple test of 
whether an isothermal total mass distribution is a good de- 
scription of other lens galaxies. In particular, we compute 
the luminosity-weighted velocity dispersion profile for a 
luminous Hernquist component in a logarithmic potential, 
which depends very mildly on the anisotropy radius, espe- 
cially at radii larger than R e . The measured velocity dis- 
persion of the only system available besides MG2016+112 
(0047+281) agrees very well with the curve computed with 
this model. As new velocity dispersion measurements of 
E/S0 lens galaxies become available, they can be readily 
compared to the expectations of this simple model to test 



8 Plus showing no evidence for rotation, since this dynamical model is inappropriate for rotating systems. 
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whether indeed the isothermal mass profile is a general 
feature of the E/SO population. 

6.2. Discussion 

In this paper and KT02 we presented the following sim- 
ple picture of the lens galaxy D in MG2016+112. The 
stellar populations are old and metal rich. A significant 
amount of dark matter is present and the mass distribution 
of the luminous and dark component are well constrained, 
reproducing all the available observations. 

However, the good agreement of the total mass profile 
with an r~ 2 power law calls for a physical explanation. 
This seems not to be a mere coincidence, given that a sim- 
ilar behavior has been observed in some local ellipticals 
with kinematic tracers at extraordinary large radii (e. g. 
Mould et al. 1990; Franx et al. 1994; Rix et al. 1997). 
Whether this galaxy formed by monolithic collapse, or 
by mergers of smaller subunits in a hierarchical picture, 
why do these processes conspire to produce an apparent 
isothermal mass distribution within R,Einst? In the case 
of adiabatic contraction, there appears to be no reason 
why gas could not continue to contract or stop earlier and 
give rise to a total mass distribution that is either steeper 
or shallower than isothermal. Thus, formation by adia- 
batic contraction seems not to be a satisfactory explana- 
tion, and it is worth seeking a more fundamental expla- 
nation. By contrast, complete violent relaxation produces 
isothermal distribution functions (Lynden-Bell 1967; Shu 
1978), and would explain quite naturally the observations. 
However, this only applies exactly for an infinite mass dis- 
tribution. When finite mass distributions are considered, 
(incomplete) violent relaxation leads to an i? 1 / 4 surface- 
brightness profile (e.g. van Albada 1982; Hjorth & Madsen 
1991, 1995), which goes as r~ 2 within the half mass radius, 
the region of interest here. For the two-component case 
(luminous and dark), the process should be such as the 
two components interact to produce the total r -2 density 
profile, while preserving the segregation of the luminous 
component in a centrally concentrated density profile (see 
e.g. the two-component models of Bertin, Saglia & Sti- 
avelli 1992; a similar outcome can be produced by merg- 



ing of galaxies with massive halos, see e. g. Barnes 1992; 
Barnes & Hernquist 1996). In the CDM scenario, it is not 
clear whether the chaotic changes in the gravitational po- 
tential during the collapse of halos satisfy the conditions 
for violent relaxation (e.g. Flores & Primack 1994 and ref- 
erences therein), although it might not be inconceivable 
that violent relaxation of the inner parts of E/SO galaxies 
is close to complete. A complete discussion of these issues 
goes beyond the aims of this paper. However, we would 
like to emphasize that whatever physical processes are in- 
volved, the mass distribution of galaxy D in MG2016+112 
seems already relaxed ~ 8 Gyr ago. 

Finally, we note that our poor knowledge of the slope of 
the mass profile of lens galaxies is the major systematic un- 
certainty in the determination of Ho from time-delays (e.g. 
Koopmans & Fassnacht 1999; Romanowsky & Kochanek 
1999; Koopmans 2001). The assumption of isothermal 
mass distribution leads to Ho ~ 60—70 kms -1 Mpc -1 (e.g. 
Koopmans & Fassnacht 1999; Koopmans 2001), in good 
agreement with other independent measurements. The re- 
sult for MG2016+112 seems to indicate that the system- 
atic error - the deviation from isothermality - might in- 
deed be very small for E/SO lens galaxies. However, the 
LSD Survey by measuring the average effective slope of 
the mass profile and its scatter, will allow us to pin down 
the largest systematic uncertainty in the global determi- 
nation of Ho from gravitational-lens time delays. 

We are grateful to E. Agol, A. Benson, G. Bertin, 
R. Blandford, L. Ciotti, R. Ellis, M. Stiavelli for the stimu- 
lating conversations and insightful comments that helped 
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